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Abstract 

We present the transverse momentum (pr) spectra for identified charged pions, 
protons and anti-protons from p+p and ci+Au collisions at ^/snn = 200 GeV. The 
spectra are measured around midrapidity (|y| < 0.5) over the range of 0.3 < p^ < 10 
GeV/c with particle identification from the ionization energy loss and its relativistic 
rise in the Time Projection Chamber and Time-of-Flight in STAR. The charged pion 
and proton+anti-proton spectra at high p^ in p+p and <i+Au collisions are in good 
agreement with a phenomenological model (EPOS) and with next-to-leading order 
perturbative quantum chromodynamic (NLO pQCD) calculations with a specific 
fragmentation scheme and factorization scale. We found that all proton, anti-proton 
and charged pion spectra in p+p collisions follow xx-scaling for the momentum 
range where particle production is dominated by hard processes (pr ~ 2 GeV/c). 
The nuclear modification factor around midrapidity is found to be greater than 
unity for charged pions and to be even larger for protons at 2 < px < 5 GeV/c. 



Key words: Particle production, perturbative quantum chromodynamics, 
fragmentation function, Cronin effect and xx-scaling. 



1 Introduction 



The study of identified hadron spectra at large transverse momentum (px) i n 
p+p collisions can be used to test the predictions from perturbative quantum 
chromodynamics (pQCD) [1]. In the framework of models based on QCD, 
the inclusive production of single hadrons is described by the convolution of 
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parton distribution functions (PDFs), parton interaction cross-sections and 
fragmentation functions (FFs). The PDF provide the probability of finding a 
parton (a quark or a gluon) in a hadron as a function of the fraction of the 
hadron's momentum carried by the parton. The FFs [2] give the probability 
for a hard scattered parton to fragment into a hadron of a given momentum 
fraction. These are not yet calculable from the first principles and hence are 
generally obtained from experimental data (e.g., e + +e~ collisions). The fac- 
torization theorem for cross-sections assumes that FFs are independent of the 
process in which they have been determined and hence represent a universal 
property of hadronization. It is therefore possible to make quantitative pre- 
dictions for other types of collision systems (e.g., p+p ). Comparisons between 
experimental data and theory can help to constrain the quark and gluon FFs 
that are critical to predictions of hadron spectra in p+p, p+A, and A+A col- 
lisions. The simultaneous study of identified hadron p T spectra in p+p and 
d+Au collisions may also provide important information on the PDFs [3] of 
the nucleus. 

The identified particle spectra in p+p and <i+Au collisions also provide refer- 
ence spectra for particle production at high px in Au+Au collisions. Moreover, 
studies of identified particle production and their ratios as a function of px in 
high-energy heavy-ion collisions have revealed many unique features in differ- 
ent px regions [4-7] and between baryons and mesons [8]. A good description 
of both identified pion and proton spectra in p+p and d+Au collisions at in- 
termediate and high px by NLO pQCD will provide a solid ground for models 
based on jet quenching [9] and quark recombination [6]. These emphasize the 
need for a systematic study of px spectra from p+p and d+Au collisions at 
the same energy as the nucleus-nucleus collisions. 

In this letter, we present the px spectra for identified pions, protons and anti- 
protons in p+p and d+Au collisions at v /snn = 200 GeV as measured by 
the STAR experiment at RHIC. The results are compared to NLO pQCD 
calculations and a phenomenological model. We also study the xx-scaling in 
p+p collisions and the nuclear modification factors in <i+Au collisions. 



2 Experiment and Analysis 

The STAR experiment consists of several detectors to measure hadronic and 
electromagnetic observables spanning a large region of the available phase 
space at RHIC. The detectors used in the present analysis are the Time Pro- 
jection Chamber (TPC), the Time-Of-Flight (TOF) detector, a set of trigger 
detectors used for obtaining the minimum bias data, and the Forward Time 
Projection Chamber for the collision centrality determination in d+Au colli- 
sions. The details of the design and other characteristics of the detectors can 
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Fig. 1. dE/dx distribution normalized by pion dE/dx at 4.5 < < 5.0 GeV/c 
and | rj | < 0.5, and shifted by ±5 for positively and negatively charged particles, 
respectively. The distributions are for minimum bias d+Au collisions. The pion, 
proton and anti-proton peak positions are indicated by arrows. 

be found in Ref. [10]. 

A total of 8.2 million minimum bias p+p collision events and 11.7 million d+Au 
collision events have been analyzed for the present study. The data set was 
collected during the years 2001 and 2003. The details of minimum bias trigger 
conditions for p+p and d+Au collisions can be found in the Refs. [11,12]. The 
minimum-bias trigger captured 95 ±3% of the 2.21 ±0.09 barn d+Au inelastic 
cross-section. The trigger efficiency was determined from a cross study of two 
sets of trigger detectors: two Zero-Degree Calorimeters (ZDCs) and two beam- 
beam counters (BBCs). The absolute cross-section is derived from a Monte 
Carlo Glauber calculation. These results are consistent with other recent mea- 
surements [13]. The trigger for the minimum bias p+p collisions required a 
coincidence measurement of the two BBCs covering 3.3 < | r] | < 5.0 [14]. 
This trigger was sensitive to color exchange hadronic and doubly-diffractive 
events; here, these are labelled "non-singly diffractive (NSD) events". Using 
PYTHIA(v6.205) [15] and HERWIG [16], it was determined that the trigger 
measured 87% of the 30.0 ± 3.5 mb NSD cross-section, which was measured 
via a vernier scan [17]. The data from TOF are used to obtain the identified 
hadron spectra for p^ < 2.5 GeV/c. The procedure for particle identification 
in TOF has been described in Ref. [18]. For p T > 2.5 GeV/c, we use data 
from the TPC. Particle identification at high pt in the TPC comes from the 
relativistic rise of the ionization energy loss (dE/dx). Details of the method 
are described in Ref. [19]. At Pt^3 GeV/c, the pion dE/dx is about 10-20% 
higher than that of kaons and protons due to the relativistic rise, resulting 
in a few standard deviations (l-3cx) separation between them. Since pions are 
the dominant component of the hadrons in p+p and d+Au collisions at RHIC, 
the prominent pion peak in the dE/dx distribution is fit with a Gaussian to 
extract the pion yield [19]. The proton yield is obtained by integrating the 
entries (Y) in the low part of the dE/dx distribution about 2.5a away from 
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the pion dE/dx peak. The integration limits were varied to check the stability 
of the results. Fig. 1 shows a typical dE/dx distribution normalized by the 
pion dE/dx at 4.5 < pr < 5.0 GeV/c and | rj | < 0.5. The Gaussian distribu- 
tion used to extract the pion yield and the pion, proton and anti-proton peak 
positions are also shown in the figure. 

The kaon contamination is estimated via either of the equations given below. 
The uncorrected proton yield is 

p = (Y-P(h-n))/(a-P) 

or 

p=(Y-pK$)/a, 

where a and (3 are the proton and kaon efficiencies from the integration de- 
scribed above, derived from the dE/dx calibration, resolution and the Bichsel 
function [19,20]. In the first case the kaon contamination is estimated through 
the yields of the inclusive hadrons (h) and pions, in case two from known yields 
from Kg measurements [19,21]. The typical values of a for a dE/dx cut slightly 
away from the proton peak position is 0.4 and the f3 values decrease from 0.2 
to 0.08 with p T in the range 2.5 < p T < 10 GeV/c. At high p T , the yields of 
other stable particles (i.e., electrons and deuterons) are at least two orders of 
magnitude lower than those of pions, and are negligible in our studies. The 
two results are consistent where STAR K$ measurements are available. The 
pr-dependence of the reconstruction efficiency, background and the systematic 
uncertainties for pions, protons and anti-protons for low in P+P and ci+Au 
collisions are described in Ref. [18]. At high p T (> 2.5 GeV/c), the efficiency 
is almost independent of in both p+p and d+Au collisions. The tracking 
efficiencies are ~ 88% and 92% in p+p and d+Au collisions, respectively. The 
difference in tracking efficiency arises because of worse vertex determination 
in p+p collisions than d+Au collisions. The background contamination to pion 
spectra for p T > 2.5 GeV/c, primarily from K$ weak decay is estimated from 
PYTHIA/HIJING simulations with full GEANT detector descriptions to be 
~ 4%. The charged pion spectra are corrected for efficiency and background 
effects. The inclusive proton and anti-proton spectra are presented with effi- 
ciency corrections and without hyperon feed-down corrections. The integrated 
A/p-ratio is estimated to be < 25% [18,21]. Additional corrections are applied 
for primary vertex reconstruction inefficiency as discussed in Refs. [11,12,18]. 
The momentum resolution is given as Ap-p/Vr = 0.01 + 0.005pT/(GeV/c) 
and has < 4% effect on the yields at the highest pr value. The spectra are 
not corrected for momentum resolution effects, but they are included in the 
systematic errors. 

The total systematic uncertainties associated with pion yields are estimated 
to be < 15%. This systematic uncertainty is dominated by the uncertainty in 
modeling the detector response in the Monte Carlo simulations. Protons from 
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Fig. 2. Midrapidity (|y| < 0.5) transverse momentum spectra for charged pions, 
proton and anti-proton in p+p and d+Au collisions for various event centrality 
classes. Minimum bias distributions are fit to Levy functions which are shown as 
dashed curves. 



hyperon (A and E) decays away from the primary vertex can be reconstructed 
as primordial protons at a slightly higher px than their true value, but with 
worse momentum resolution. This results in an uncertainty of the inclusive 
proton yield of ~ 2% at p T = 3 GeV/c and ~ 10% at p T = 10 GeV/c. For 
proton and anti-proton yields at high px an additional systematic error arises 
from the uncertainties in the determination of the efficiencies, a and /3, under 
a specific dE/dx selection for integration. This is due to the uncertainties 
in the mean dE/dx positions for protons and kaons. The total systematic 
uncertainty in obtaining the proton and anti-proton yields for px > 2.5 GeV/c 
increases with p T from 12% to 23% (at pr = 10 GeV/ c) in both p+p and d+Au 
collisions. The errors shown in the figures are statistical, and the systematic 
errors are plotted as shaded bands. In addition, there are overall normalization 
uncertainties from trigger and luminosity in p+p and d+Au collisions of 14% 
and 10%, respectively [11]. These errors are not shown. 

Figure 2 shows the invariant yields of charged pions, protons and anti-protons 
for the pt range of 0.3 < px < 10 GeV/c in minimum bias p+p collisions 
and for various centrality classes in d+Au collisions. The yields span over 
eight orders of magnitude. The minimum bias distributions are fit with a 
L evy distr ibution [22] of the form 2 J^ prdy = (1+(mr _t o)/nTr , where m T = 

Jp\ + tuq and mo is the mass of the hadron. The Levy distribution essentially 
takes a power-law form at higher px and has an exponential form at low pt- 
For the p and p spectra, fit with a power-law function gives a worse x 2 /ndf 
compared to the fit with the Levy function. For d+Au collisions the x 1 jndj for 
the power-law fit to p(p) spectra is 68.55/20(86.77/20) and the corresponding 
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Fig. 3. Nuclear modification factors, RdAu, for charged pions ir + + tt~ and p+p 
at \y\ < 0.5 in minimum bias and 0-20% central d+Au collisions. For comparison 
results on inclusive charged hadrons (STAR) from Ref. [11] at \r]\ < 0.5 are shown 
by dashed curves. The first two shaded bands around 1 correspond to the error due 
to uncertainties in estimating the number of binary collisions in minimum bias and 
0-20% central d+Au collisions respectively. The last shaded band is the normaliza- 
tion uncertainty from trigger and luminosity in p+p and d+Au collisions. 

value for the fit with the Levy function is 21.19/20(26.4/20). 



3 Nuclear modification factor 

The nuclear modification factor (-RdAu) can be used to study the effects of cold 
nuclear matter on particle production. It is denned as a ratio of the invariant 
yields of the produced particles in d+Au collisions to those in p+p collisions 
scaled by the underlying number of nucleon-nucleon binary collisions. 



where (Nun) is the average number of binary nucleon-nucleon (NN) collisions 
per event, and (A^bin)/o"pp el is the nuclear overlap function TA_(b) [11,12]. The 
value of up" 61 is taken to be 42 mb. 

The left panel of Fig. 3 shows -RdAu values for charged pions ((7r + +7r~)/2) in 
minimum bias and 0-20% central collisions at |y| < 0.5. The -RdAu for 0-20% 
central collisions are higher than -RdAu for minimum bias collisions. The result 
-RdAu > 1 indicates a slight enhancement of high p T charged pion yields in 
d+Au collisions compared to binary collision scaled charged pion yields in 
p+p collisions within the measured (y, px) range. The right panel of Fig. 3 
shows the -RdAu of baryons (p+p) for the minimum bias collisions at |y| < 0.5. 
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Fig. 4. Ratio of tt~ /ir + , p/p, p/ir + , p/ir~ at midrapidity (\y\ < 0.5) as a function of 
PT in p+p minimum bias collisions. For comparison the results from lower energies 
at ISR [26] and FNAL [27] are also shown for p/ir + and p/it~ ratios. The dotted 
curves are the results from PYTHIA. The shaded bands below the ir~ /tt + and p/p 
ratios are the point-to-point correlated errors in the yields associated with the ratio. 

The -RdAu for p+p is again greater than unity for p^ > 1.0 GeV/c and is larger 
than -RdAu for charged pions. The -RdAu of pions for 2< p^ < 5 GeV/c is 1.24 
± 0.13 and that for p+p is 1.49 ± 0.17 in minimum bias collisions. Identified 
hadron -RdAu are sensitive to nuclear modification of the PDF from processes 
such as nuclear shadowing and parton saturation as well as to transverse 
momentum broadening, energy loss in cold nuclear matter and hadronization 
through recombination, thereby further constraining the models [23]. 



4 Particle ratios 



The particle ratios at midrapidity as a function of pi for p+p and <i+Au 
minimum bias collisions are shown in Figs. 4 and 5 respectively. Correlated 
errors are shown as the shaded bands below the data points. The 7r~/7r + -ratio 
has a value ~ 1 and is independent of p^ in both p+p and <i+Au collisions. 
The p/p-ratio for p+p collisions is also independent of px within the range 
studied and has a value of 0.81 ± 0.1 at 2.5 < px < 6.5 GeV/c. However, in 
d+Au collisions we observe a clear decrease of p/p for p T >6 GeV/c. In quark 
fragmentation, the leading hadron is more likely to be a particle rather than 
an anti-particle, and there is no such preference from a gluon jet. A decrease 
in the antiparticle/particle ratio with p^ would then indicate a significant 
quark jet contribution to the baryon production. It is, however, not clear 
whether the same effect exists in p+p collisions or whether the decrease of 
p/p is due to additional nuclear effects in d+Au collisions. Calculations from 
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Fig. 5. Same as Fig. 4 for d+Au minimum bias collisions. For comparison the 
p/7r + -ratio from lower energies at FNAL [27] are shown. 

PYTHIA(v6.319) predict somewhat more prominent pi-dependence [15]. 

At RHIC, the p/n + and p/ir~ ratios increase with p T up to 2 GeV/c and 
then start to decrease for higher p^ in both p+p and o?+Au collisions. The 
p/7r~-ratio rapidly approaches a value of 0.2, which is between the values 
in e + +e" collisions for quark and gluon jets [24,25]. The p/n + and p/n~ 
ratios from PYTHIA are constant at high p T in contrast to a decreasing trend 
oberserved in the data. The p/7r + -ratios in p+p collisions compare well with 
results from lower energy ISR and FNAL fixed target experiments [26,27]. 
Meanwhile, p/7r~-ratios at high p^< have a strong energy dependence with 
larger values at higher beam energies. In <i+Au collisions the p/7r + -ratio at 
high p T is lower for p+A collisions at FNAL energy than at RHIC. 



5 Comparison to NLO pQCD and model calculations 



In Fig. 6 we compare (n + + 7r~)/2 and (p+p)/ 2 yields in minimum bias 
p+p and d+Au collisions at midrapidity for high p T to those from NLO 
pQCD calculations and the phenomenological parton model (EPOS) [28]. 
The results from EPOS agree fairly well with our data for charged pions 
and proton+anti-proton in p+p and d+Au collisions. The NLO pQCD results 
are based on calculations performed with two sets of FFs, the Kniehl- Kramer- 
Potter (KKP) [29] and the Albino- Kniehl- Kramer (AKK) set of functions [30]. 
The factorization scale for all the NLO pQCD calculations shown is for fi = 
Pt- The charged pion data for pi > 2 GeV/c in p+p collisions are reasonably 
well described by the NLO pQCD calculations using the KKP and AKK set of 
FFs. A similar observation for n°s using KKP FFs was made by the PHENIX 
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Fig. 6. Midrapidity invariant yields for (ir + + ix )/2 and {p+p)/2 at high p^ for 
minimum bias p+p and d+Au collisions compared to results from NLO pQCD 
calculations using KKP [29] (PDF: CTEQ6.0) and AKK [30] (PDF: CTEQ6M) sets 
of fragmentation functions and results from the EPOS model [28]. The PDFs for 
d and Au-nucleus are taken from Refs. [31] and [32] respectively. All results from 
NLO pQCD calculations are with factorization scale is /x = p^. 



Collaboration [33]. For d+Au collisions NLO pQCD calculations with KKP 
FFs are consistent with the data for p^ > 4 GeV/c while those with AKK 
FFs underpredict the measured charged pion yields. 

The proton+anti-proton yield at high p^ in p+p and d+Au collisions is much 
higher than the results from NLO pQCD calculations using the KKP set 
of FFs and lower compared to calculations using AKK FFs. The relatively 
better agreement of NLO pQCD calculations with AKK FFs compared to 
those with KKP FFs for proton+anti-proton yields shows the importance of 
the flavor-specific measurements in e + +e~ collisions in determining the FFs 
for baryons. One may further improve the NLO pQCD calculations by an 
all-order resummation of large logarithmic corrections to the partonic cross- 
sections [34]. 



6 Scaling of particle production 



The invariant cross- sect ions of inclusive pion production in high energy p+p 
collisions have been found to follow the scaling laws [36] : 

rf 3 a 1 „ , . ^d 3 a 1 
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Fig. 7. XT-scaling of pions, protons and anti-protons. The data from other exper- 
iments are from the following references, FNAL : Refs [27,35], ISR : Ref. [26], 
PHENIX : Ref. [33], and UA2 [40]. The inset shows the mx-scaling of the invariant 
yields for charged pions and protons+anti-protons in p+p and d+Au collisions. 

where xt = 2p^/y/s and / (xt) and g (xt) are some functions of x<r- Similar 
scaling has been observed in e + +e~ collisions, but without the a/s™ or p T fac- 
tor [37]. The value of the power n ranges from 4 to 8 [38]. In the general scaling 
form ~ 1/pt> n depends on the quantum exchanged in the hard scattering. In 
parton models, it is related to the number of point-like constituents taking an 
active role in the interaction. The value reaches 8 in the case of a quark-meson 
scattering by exchanging a quark. With the inclusion of QCD, the scaling law 
follows as ~ l/i/i™, where n becomes a function of xt and y/s. The value of n 
depends on the evolution of the structure function and FFs. n=4 is expected 
in more basic scattering processes (as in QED) [38,39]. 

Figure 7 shows the xi-scaling of pions, protons and anti-protons. The value of 
n obtained for the scaling with yfs 1 of the invariant cross-section is 6.5 ± 0.8. 
The STAR data covers the range 0.003 < x-r < 0.1. The data points deviate 
from the scaling behavior for p^ < 2 GeV/ c for pions and protons, which could 
be interpreted as a transition region from soft to hard processes in the particle 
production. The deviations start at a higher p^ for the anti-protons. The 
available data on pion and proton invariant cross-sections at various center- 
of-mass energies [26,27,33,35,36,40] for p T > 2 GeV/c are compiled and fitted 
using the function 4r (1 — x^) m . The value of n ranges from 6.0 to 7.3 for 
v /snn between 19 GeV and 540 GeV, while that for m ranges between 13 
and 22. The average value of n for pions is 6.8 ± 0.5 and that for protons 
and anti-protons is 6.5 ± 1.0. The variations in n and m values may lead 
to differences in details of scaling behaviour at different energies when the 
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cross-section is multiplied by 1/pt [41]- This feature is not observed in the 
scaling shown in Fig. 7 due to the data spanning several orders of magnitude. 
The inset of Fig. 7 shows the m T scaling at px < 2 GeV/c, consistent with 
possible transition between soft and hard processes at around p T — 2 GeV/c. 
The mT-scaling also indicates that flow effects in p+p and d+Au collisions are 
negligible [4,5]. The presented data suggests that the transition region from 
soft to hard physics occurs around px~2 GeV/c in p+p collisions. 



7 Summary 

We have presented transverse momentum spectra for identified charged pi- 
ons, protons and anti-protons from p+p and d+Au collisions at v /snn = 200 
GeV. The transverse momentum spectra are measured around midrapidity 
(|y| < 0.5) over the range of 0.3 < pi < 10 GeV/c with particle identification 
from the ionization energy loss and its relativistic rise in the Time Projection 
Chamber, as well as the Time-of-Flight in STAR. The following conclusions 
can be drawn from the present study: (a) The nuclear modification factor 
around midrapidity is enhanced in d+Au collisions to about 1.5 for pions and 
to about 2 for protons and antiprotons at intermediate pr (2 < Pt < 5 GeV/c). 
(b) Identified particle ratios were measured up to pt of 7 GeV/c in p+p and 
10 GeV/c in d+Au reactions. Their dependence on species, px and collisions 
energy was shown to be sensitive to the relative contributions from quark and 
gluon fragmentation as well as to their fragmentation functions, (c) The NLO 
pQCD calculations describe the high p T data for charged pions reasonably well 
in p+p collisions and d+Au collisions. In general, baryon production has his- 
torically been difficult to describe by pQCD and hadronization [39,42]. Use of 
the recently published AKK FFs results in a much improved description of the 
measured p and p spectra, (d) The proton and pion spectra in p+p collisions 
follow ^T-scaling with a beam-energy dependent factor ~ ^snn 6 ' 5 above px~ 
2 GeV/c. The pion and proton spectra follow transverse mass scaling for 
< 2 GeV/c 2 in both p+p and d+Au collisions, suggesting the transition region 
from soft to hard process domination occurs at pt~ 2 GeV/c in these collision 
systems. The measurements presented in this paper provide better constraints 
on jet quenching and quark recombination models which are presently the best 
candidates for explaining particle production in the intermediate px region. 
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